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A linear stability analysis of metallic nanowires is performed in the free-electron 
model using quantum chaos techniques. It is found that the classical instability of 
a long wire under surface tension can be completely suppressed by electronic shell ef- 
fects, leading to stable cylindrical configurations whose electrical conductance is a magic 
number 1, 3, 5, 6,... times the quantum of conductance. Our results are quantitatively 
consistent with recent experiments with alkali metal nanowires. 



I. INTRODUCTION 

Plateau's celebrated study of the stability of bodies under the influence of surface tension estab- 
lished a fundamental result of classical continuum mechanics: a cylinder longer than its circumference is 
energetically unstable to breakup. Here we consider a quantum mechanical generalization of Plateau's 
study — the stability of a metallic nanowire. Methods from the study of quantum chaos ||,Q are used 
to examine the energetics of a free electron gas in a cylindrical jellium filament, and show that an 
infinite filament can be stable against all axisymmetric perturbations if its electrical conductance G 
[in units of Go = 2e 2 /h = (12, 906.4 Ohm) -1 ] belongs to a set of "magic numbers" n = 1,3,5,6, .. . 
and is otherwise unstable. Here e is the charge of an electron and h is Planck's constant. Such magic 
numbers are analogous to those found in the shell model of atomic nuclei Q and in metal clusters 
[^,^) . Our stability analysis elucidates and confirms the stability properties of simple metal nanowires 
conjectured by Yanson, Yanson, and van Ruitenbeek |j] based on their observation of shell structure in 
the conductance statistics of sodium nanowires. 

Structural multistability of metallic nanowires was previously postulated based on classical molecular 
dynamics simulations However, the magic numbers observed in conductance histograms in alkali 

metals [p|]To|] are clearly an electronic shell effect, as shown below, and cannot be explained with classical 
molecular dynamics. The common occurence of multistability in two such radically different models 
likely stems from the fact that both models introduce an additional length scale (the Fermi wavelength 
\p in the free electron model, the atomic diameter in molecular dynamics simulations), which leads to 
commensurability effects. 

The properties of nanowires formed from simple monovalent metals, like the bulk properties of these 
materials nj, are determined largely by the delocalized conduction electrons. A free electron jellium 
model, treating the electrons as a noninteracting Fermi gas confined within the wire by hard-wall 
boundary conditions, provides an intuitive and even quantitative explanation of observed quantities 
like conductance jl2]-[l4| , force |l3|~|l7|] , and shot noise JlJ] . To examine energetic stability, we consider 
a weakly deformed cylinder, and find its thermodynamic potential to quadratic order in the amplitude 
of the deformation. A contribution proportional to the filament area, naturally identified as the surface 
tension, is only one of several competing terms in the free energy, which consists of terms which vary 
smoothly with the geometry of the filament and an oscillatory contribution directly connected with 
the discrete energy levels that are solutions to Schrodinger's equation. The smooth terms appear in 
descending powers of length (proportional to volume, surface area, mean curvature, etc.), and are quite 
analogous to those found in the study of classical wave equations in curved domains Jl9| and the related 
problem of classically screened Coulomb (Dcbye-Hiickel) interactions of curved surfaces in an electrolyte 
p0| |. The oscillatory part of the free energy, in particular, alters Plateau's classical stability analysis in 
an essential way. The present work is thus fundamentally distinct from that which has dealt with the 
quantum mechanical origin of surface tension itself in metallic fluids pl[ , as well as those which consider 
quantum mechanical corrections to classical surface tension due to the quantization of capillary waves 

Long gold nanowires suspended between gold electrodes have been produced and imaged with a 
transmission electron microscope by Kondo and Takayanagi p^j2^ ]; in particular, they have observed 
wires which are stable and almost perfectly cylindrical. The wires of Ref. p3j23] are of finite length, 
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FIG. 1. Deformation of wavevector q of a cylindrical nanowire. 



and attached to electrodes at either end. It seems reasonable to consider an infinite wire as a theoret- 
ical starting point, and to assume that the length of the wire acts as a cut-off: deformations with a 
wavelength longer than the wire length cannot occur. Thus one does not have to model the boundary 
conditions at the ends of the wire, which would introduce additional parameters. 

To be specific, we examine an infinite cylindrical wire and its sole classically unstable deformation — an 
axially symmetric one (see Fig. []]). Any such deformation can be written as a Fourier series: 



R{ Z )=R \l + J dqb(q)cos(qz + (p(q))\, (1) 

where R(z) is the radius at position z, Ro is the unperturbed radius, b(q) is the (infinitesimal) pertur- 
bation coefficient, and <fi(q) is an arbitrary phase shift. The coefficients b{q) are chosen such that the 
total volume of the wire is unchanged by the deformation. Other physically reasonable constraints jl7| 
are also possible, but lead to similar results. 

The metallic nanowire is an open system, connected to macroscopic metallic electrodes at each end 
J|[l0]j2^ |2^] . Therefore the change of the grand canonical potential fl under the perturbation determines 
its stability, fl is related to the electronic density of states g(E) by 

n = -k B T^ dEg{E)\n(l + e-( E -ti/ kBT } , (2) 

where k B is Boltzmann's constant, T is the temperature, and [i is the electronic chemical potential 
specified by the macroscopic electrodes. Our aim is to expand Q up to second order in the coefficients 
b(q) characterizing the deformation. As we will show, this yields 



fi[6] = fi[0] + / dqa(q)[b(q)} 2 + 0(b 3 ) . (3) 
Jo 

The change in the grand canonical potential is of second order in b and contributions from deformations 
with different q decouple. If the prefactor a(q) is negative for any value of q, then fl decreases under 
the deformation and the wire is unstable. 

For an infinite cylindrical wire, the transverse motion is quantized, leading to the formation of discrete 
electronic subbands. The total density of states is the sum of the contributions from each subband (see 
Fig. |^): every subband begins to contribute at a threshold energy equal to the energy of its quantized 
transverse motion with a sharp spike, falling off smoothly for increasing energy. If the Fermi energy Ep 
lies near one of these sharp peaks, certain small deformations of the wire can dramatically increase the 
density of states. According to (||), this lowers the grand canonical potential, leading to an instability. 
On the other hand, if there is no subband threshold sufficiently close to Ep, the density of states will 
instead decrease with any deformation, implying the existence of stable regions intervening between the 
instabilities associated with the opening of each subband. 



II. QUANTUM CHAOS APPROACH 

In order to examine this picture quantitatively, we use a semiclassical approach which enables 
g(E) to be split into a smooth average contribution g(E), referred to as the Weyl contribution, and an 
oscillatory part 8g{E), whose average value is zero: 
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FIG. 2. Density of states g(Ep) of a cylindrical wire (upper diagram) and stability coefficient a (lower 
diagram) versus the radius Rq of the unperturbed wire. The wavevector of the perturbation is qRo = 1. Dashed 
curve: Weyl contribution to a. 



g(E)=g(E)+Sg(E). (4) 

The Weyl term g contains terms proportional to the volume of the nanowire, and to the area and 
curvature of its surface: g(E) = (1/ 'Ep)g(E / 'Ep), where 

2.1/2 i -1/2 

and kp = 2tt/\p is the Fermi wavevector. The volume V, surface area S, and integrated mean curvature 
K can be calculated for arbitrary perturbations by simple geometric considerations. The oscillatory 
part Sg of the density of states is a quantum correction, and may be calculated in the semiclassical 
approximation as a sum over all periodic classical orbits of the system P,^, p^| , p7| , p9| -|3~H . 

Let us first consider an undcformcd cylindrical wire. Each periodic orbit lies in a plane transverse to 
the wire's axis, and there is a correspondence to the orbits Q in a circular billiard (see Fig. ||). Note 
that the action S vw of each orbit is invariant under both translations of the orbit parallel to the z-axis 
and rotations about the z-axis. In a system with continuous symmetries, when taking the trace of the 
electron Green's function to calculate the density of states, one must first integrate exactly over these 
symmetries [| 29 before employing the semiclassical (stationary phase) approximation. For electrons 



in a cylinder of length L, one obtains the following trace formula 

tiiL \ \ fvwL vw ^ f S VW (E) 3w7r 



w— 1 v— 2w x ' 

where m is the electron mass, v and w are defined in Fig. |[ f vw = 1 + 6{v — 2w) counts the discrete 
degeneracy of the orbit under time reversal, L vw = 2vRq sm^irw/v) is the length of a periodic orbit, 
and the action S vw /h = kpL vw ^J E / Ep . 

When the wire is deformed, the translational symmetry is broken, and we use a semiclassical pertur- 
bation theory [ |30|]3l"| ] to calculate 5g. This theory assumes that for small perturbations, the amplitudes 
and topology of the unperturbed orbits can still be used in but that the lengths of the orbits — and 
hence their actions — change. The trace formula becomes 

§ 9( E ) = ^EE ^ Re e *(W*-3W2) / dz {Lvw + A£ w)e <AW» , (7) 

^ iu — 1 v — 2w L ^ J 
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FIG. 3. Classical periodic orbits of an electron in a plane perpendicular to the z-axis, labeled (v,w), where v 
is the number of vertices and w the winding number. 



where AS vw /h — kpAL vw ^E/E F , and 

AL VW = 2vsm(TTw/v)(R(z) - R ), (8) 

which may be expressed in terms of the perturbation coefficients using Eq. (Q). Combining this result for 
Sg with (g), it is straightforward to calculate the density of states up to second order in the coefficients 
b(q) for a deformed wire. The result can be integrated to obtain the grand canonical potential, which 
indeed has the expansion (ph. 



III. STABILITY ANALYSIS 



Let us first discuss the stability of a nanowire at zero temperature. Fig. |2| shows the stability coefficient 
a (lower diagram) and the density of states at the Fermi energy (upper diagram) as functions of the 
radius of the unperturbed wire. The wavelength of the perturbation was taken to be qRq = 1, the critical 
wavelength for stability in Plateau's classical analysis of a body under the influence of surface tension. 
In addition to surface tension, the present model for f2 has a curvature energy, which enhances the 
instability for small Rq, and an oscillatory component associated with the opening of successive subbands 
as Rq increases. As discussed above, a has sharp negative peaks — indicating strong instabilities — at the 
subband thresholds, where the density of states is sharply peaked. Under surface tension and curvature 
energy alone (dashed curve in Fig. ||), the wire would be slightly unstable at the critical wavevector 
qRo = 1, since the curvature term is negative. However, the quantum correction is positive in the 
regions between the thresholds to open new subbands, thus stabilizing the wire. Since the oscillatory 
contribution to a is independent of q, we find that regions of stability persist for arbitrarilly long 
wavelength perturbations, indicating that an infinitely long cylindrical wire is a true metastable state if 
the radius lies in one of the windows of stability. 

With these results, we can construct the zero temperature stability diagram for the wire [see Fig. 
0(a)]. In contrast to Plateau's classical stability analysis, an additional quantum length scale arises 
here, namely the Fermi wavelength Xp- The stability problem is now determined by two dimcnsionlcss 
parameters: qRo and fcpi?o- In Fig- 0j regions of instability, where the coefficient a(q) is negative, 
are shaded grey, while the stable regions are shown in white. Note that many of the white regions of 
stability persist all the way down to q — 0. The multistability of the system, indicated by the alternating 
stable and unstable stripes, reflects commensurability effects between Xp and Rq. 

We remark that in addition to the axially symmetric modes considered here, which are the sole 
unstable modes in the classical limit, perturbations which break axial symmetry may also become 
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FIG. 4. Stability diagram for cylindrical nanowires at two different temperatures. White areas are stable, 
grey unstable to small perturbations. The quantized electrical conductance values G of the stable configurations 
are indicated by bold numerals in (a), with Go = 2e 2 /h. Right vertical axis: corrected Sharvin conductance 
Gs- Dotted curve: stability criterion in the Weyl approximation. 

unstable near the subband thresholds. However, these Jahn- Teller like modes p5[ are likely less unstable 
than the axisymmetric modes, and will not destroy the regions of stability shown in Fig. ^(a). 

IV. CONDUCTANCE MAGIC NUMBERS 

The electrical conductance G of a perfect cylindrical nanowire is quantized p2| in units of Go = 2e 2 /h. 
The quantized conductance values of the stable cylindrical configurations are indicated by bold numerals 
in Fig. ^(a). For comparison, the right vertical axis of the figure is labeled with the corrected Sharvin 
conductance lna| 



G s = 



k F R 



1 - 



(9) 



which gives a smooth approximation to G/Gq. The conductance values of the stable configurations 
are somewhat analogous to the magic numbers of enhanced stability in atomic nuclei [Q and metal 
clusters QJ^]. An important distinction is that the magic numbers in nuclei and clusters refer to the 
number of fermions in a finite system, while we consider an infinite, open system, with magic numbers 
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describing the number of conducting transverse modes |13| which hold the wire together. The number 
of conducting modes is approximately equal to the number of atoms which fit within a cross section of 
the wire. 

The integer magic numbers identified here must be distinguished from the concept of "magic wire 
configurations" proposed by van Ruitenbeek and others [|l5 32 3j|. The latter represented discrete 
minima of the energy per unit volume of a jellium wire as a function of its cross-sectional area, but the 
stability with respect to deformations with q > was not considered. A similar scenario was advanced 
in Ref. based on a partial summation of the periodic orbits in the expression for itself. While 
the theoretical results of Refs. ||,[l5|,[32],[33| seem to imply that stability occurs only for a discrete set 
of radii, of measure zero, our analysis finds stability with respect to small perturbations over broad 
intervals of radius. Furthermore, the very existence of discrete energetic minima of the type discussed 
in Refs. |(||l5|,[32|,[53| depends sensitively on the numerical value of the surface tension in the model; 
for example, they do not occur whatsoever in the free-electron model at constant volume for G > Gq. 
In contrast, the finite windows of stability in our analysis are robust with respect to variations of the 
surface tension. 



V. COMPARISON TO EXPERIMENTS IN ALKALI METALS 



The sequence of magic numbers G/Gq — 1, 3, 5, 6, ... is consistent with the observation of conductance 
quantization in alkali metal point contacts [jlO| . Recently, conductance histograms for sodium nanowires 
with pronounced peaks up to G/Gq ~ 100 were obtained by Yanson et al. Ejj. They argued that 
these peaks could not be understood based on conductance quantization alone, but rather reflected 
energetically preferred wire configurations. In order to construct a theoretical conductance histogram, 
we need to know the a priori probability of occurence of a nanowire of a given cross-sectional area. 
Here, we make the simplest hypothesis: that nanowires of different cross- sectional areas occur with 
equal a priori probability if stable, and with zero probability otherwise. On this hypothesis, the relative 
probability of observing a contact with a given quantized conductance value is proportional to the width 
AGs of the corresponding stable region, Gs being a dimensionless measure of the contact area. 

The conductance histogram from Rcf. || taken at a temperature of T = 80K is reproduced in our Fig. 
||(a). For comparison, the theoretical magic numbers at the same temperature are plotted as vertical 
bars, with height equal to the width AGs of the corresponding stable region. Unlike the idealized 
wires in our analysis, the experimental wires are of finite length, and contain imperfections. Thus 
the peaks in the experimental histogram are shifted |l4| below the theoretical integer values due to 
backscattering, and are broadened |l4| due to tunneling, disorder-induced conductance fluctuations, 
and inelastic processes. Furthermore, the relative heights of the peaks may be influenced by dynamical 
as well as energetic effects; in particular, the peak at G = Go, which is quite pronounced at the lowest 
temperatures IgJlQl, is absent from the experimental histogram at T = 80K — presumably indicating 
that thermally activated processes lead to the rupture of this metastable configuration. 

While the first four theoretical peaks at G/Gq — 1,3,5,6 can be identified unambiguously with the 
narrower peaks in the low temperature data of Refs. [pUlQ], it is not entirely clear whether one can match 
up the broad peaks in the 80K experimental histogram with individual theoretical magic numbers. In 
particular, there appears to be fine structure in the theoretical histogram which is not present in the 80K 
histogram, possibly because it is obscured due to broadening. In order to see whether the theoretical 
histogram nonetheless correctly describes the gross features of the experimental histogram, it is useful 
to take the Fourier transforms of the two histograms, singling out the contributions of the shortest 
periodic orbits. 

The actions of the periodic orbits are proportional to the radius of the wire, and hence are approxi- 
mately proportional to the square root of the conductance, by Eq. (^J). The contributions of the various 
periodic orbits to the conductance histogram can thus be extracted |3~i| ] by taking a Fourier transform 
with respect to y/ G/Gq. The Fourier power spectrum of the theoretical conductance histogram is shown 
as a solid curve in Fig. ||(b), where all magic numbers up to G/Gq = 207 were included. It shows clear 
peaks corresponding to the three shortest periodic orbits. For comparison, the experimental Fourier 
spectrum for sodium nanowires at 90K obtained by Yanson et al. ]34| is shown as a dashed curve. The 
experimental spectrum contains two broad peaks, one centered at the frequency of the diameter orbit, 
and a second which spans the frequencies of the triangle and square orbits. The experimental peaks 
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FIG. 5. Magic numbers: theory and experiment, (a) Solid curve: conductance histogram for sodium 
nanowires at T = 80K = 0.002Tf, reproduced from Ref. M. Vertical bars: quantized conductance values of the 
metastable nanowire configurations in the free electron model at the same temperature. The height of each bar 
is equal to the width AGs of the corresponding stable region [c.f. Fig. ^(a)]. (b) Solid curve: Fourier power 
spectrum of the theoretical conductance histogram [vertical bars in (a)], displaying the dominant contributions 
of the three shortest periodic orbits in Eq. (Q). Dashed curve: Fourier power spectrum of the conductance 
histogram for sodium nanowires at T = 90K, reproduced from Ref. 

Note that the theoretical Fourier 

spectrum does change significantly between 80K and 90K. 



in Fig. g](b) are broader than the theoretical peaks since the oscillatory structure in the experimental 
histogram is damped for G/Gq > 50. The widths of the theoretical peaks, in contrast, are determined 
by the maximum conductance included, max{G/Go} = 207; increasing this cutoff decreases the widths 
of the peaks. We have simply taken a cutoff larger than the maximum measured conductance, so that 
the cutoff does not introduce any artifacts in the theoretical Fourier spectrum. The overall vertical 
scale of the power spectra shown is arbitrary, so a comparison between theory and experiment must be 
based on the relative spectral weights of the various peaks. In both spectra, the weight of the diameter 
orbit is roughly half the combined spectral weights of the triangle and square orbits — a rather good 
agreement between theory and experiment. 



VI. CLASSICAL LIMIT 

At zero temperature, the pattern of stable regions separated by unstable stripes shown in Fig. ||(a) 
continues up to arbitrarilly large radii. However, at any finite temperature T, the quantum oscillations 
in a are smoothed out, and the classical stability criterion qRq > 1 is recovered asymptotically for 
sufficiently large radii. The crossover from the T = result to the classical limit occurs when fcgT ~ 
Ep(Go/G), i.e., when the thermal energy fcgT is comparable to the average transverse level spacing. 
Fig. [|(b) shows the stability diagram for T/Tp — .05, where Tp = Ep/ks is the Fermi temperature. 
One sees that the stability boundary indeed begins to cross over to the classical line qRo = 1 for 
G s > 20. 

In Fig. f|(b), there are no true metastable configurations with Gs > 25 ~ Tp/T, indicating that 
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all thicker wires would be dynamically unstable (like a column of fluid) at this temperature, once the 
electronic shell effects have been smoothed out. However, Tp = 3.75 x 1 4 K in sodium, so multistability 
from electronic shell effects can be expected to occur in sodium contacts with G/Gq < 125 up to at 
least 300K. 



VII. DISCUSSION 



It should be pointed out that thermal averaging is not the only mechanism which can suppress 
electronic shell effects. Disorder also tends to smooth out the sharply peaked structure in the density 
of states, so that one can expect a reduction of shell effects when the diameter of the wire exceeds 
the mean free path. Furthermore, the tendency of the positive ions to order themselves into regular 
arrays [^3| , p4| will certainly affect the stability of metallic nanowires. Indeed, pioneering theoretical 
investigations [@-[| of the dynamics of nanowires focused exclusively on the arrangement of the ions. 
Based on the relative importance of electronic shells and crystal structure in metal clusters ||, one 
would expect electronic shell effects to dominate the energetics of very thin wires, particularly in the 
alkali metals, with crystal structure becoming increasingly important for thicker wires, and for metals 
where the bonding is more directional. 

Although conductance histograms of gold nanowires [p6| do not exhibit the sequence of magic numbers 
predicted here and observed in alkali metals, there is some evidence that gold nanowires can otherwise be 
adequately described using the free electron model Jl3| , ^4|Jl^ ] . It is thus worthwhile to speculate about 
a possible electronic origin for the remarkable stability of wires of individual gold atoms. Linear chains 
composed of four to seven gold atoms suspended between two gold electrodes, with a conductance 
G = Go, were found to be stable in the laboratory for hours at a time p7ps| ]. Given that such a 
configuration has an enormous surface energy, its stability is at first sight surprising. However, in our 
free electron model, we find that an infinitely long wire with a conductance of Go is indeed stable with 
respect to small perturbations. 

Our results may also be relevant to explain the observed stability of gold wires with larger radii 
|^3|,^4). The atomic-scale structure of these wires exhibits a remarkable diversity, ranging from helical, 
multishell structures |24| to crystalline structures with surface reconstruction [^3| ; but a common factor 
in the observed stable structures is that they are almost perfectly cylindrical. A direct comparison to 
our theoretical stability analysis would be facilitated by measurements of the conductance of the stable 
wires, which have not yet been carried out. However, a rule of thumb is that the conductance G/Go of 
a metal wire made of monovalent atoms is roughly equal to the number of atoms in the cross section. 
Thus the hexagonal prism structure with a cross section of 30 atoms determined in Ref. [^3| does 
correspond to a stable configuration in our analysis [vertical bar at G/Gq = 30 in Fig. ||(a)]. Similarly, 
the thinnest of the helical wire structures determined in Ref. |m| , with a cross section of 8 atoms, also 
corresponds to a stable structure in our analysis, while the two thickest wire structures determined in 
Ref. PJ], with cross sections of 22 and 24 atoms, respectively, appear to straddle the predicted island of 
stability at G/Gq, = 23. Further work is needed to elucidate the interplay of atomistic stacking effects 
and electronic shell effects in these structures. 

Finally, let us comment on the dynamical evolution of a nanowire under elongation or compression. 
Consider stretching a nanowire that is initially in a stable configuration (white areas in Fig. [§. Under 
elongation, the radius of the wire decreases, so that one moves downward on the stability diagram. 
When a stability boundary is encountered, it becomes energetically (and dynamically) favorable for 
the wire to deform spontaneously, until another stable configuration of smaller radius is reached, thus 
causing the conductance to jump abruptly from one magic number to a smaller one, and conversely 
under compression. This scenario is consistent with the claim p5[ that the structure of a metallic 
nanowire undergoes a sequence of abrupt changes as a function of elongation or compression. The 
finite widths of the unstable tongues in Fig. [| also provides a possible explanation for the hysteresis 
p5] | observed in the conductance as a function of elongation: the critical radius at which the wire's 
conductance jumps between neighboring magic numbers is different, depending on whether the tongue 
is approached from above or below, i.e., depending on whether the wire is stretched or compressed. 
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